Magnetic anisotropy of the spin ice compound Dy 2 Ti 2 07 
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We report magnetization and ac susceptibility of single crystals of the spin ice compound Dy2Ti207. 
Saturated moments at 1.8 K along the charasteristic axes [100] and [110] agree with the expected 
values for an effective ferromagnetic nearest-neighbor fsing pyrochlore with local (111) anisotropy, 
where each magnetic moment is constrained to obey the 'ice-rule'. At high enough magnetic fields 
along the [111] axis, the saturated moment exhibits a beaking of the ice-rule; it agrees with the 
value expected for a three-in one-out spin configuration. Assuming the realistic magnetic interaction 
between Dy ions given by the dipolar spin ice model, we completely reproduce the results at 2 K by 
Monte Carlo calculations. However, down to at least 60 mK, we have not found any experimental 
evidence of the long-range magnetic ordering predicted by this model to occur at around 180 mK. 
Instead, we confirm the spin freezing of the system below 0.5 K. 

PACS numbers: 75.40. Cx, 75.10.Hk, 05.50.+q, 75.40.Mg 



I. INTRODUCTION 

The observed absence of magnetic ordering in the Ising 
pyrochlore magnet Ho2Ti207 had a strong impact on the 
present study of frustrated systems!. This compound 
was known to have a ferromagnetic (FM) Curie- Weiss 
temperature $cw ~ 1-9 K, and magnetic Ho 3+ ions which 
possess a strong Ising anisotropy along the local (111) 
axis. Previous to this discovery, it was expected that 
absence of ordering would take place in the pyrochlore 
lattice only with antiferromagnetic Heisenberg spinsBB. 

Harris et al. showed the possibility of a macroscopi- 
cally degenerate ground state in Ho2Ti207 by explaining 
the multiplicity of the possible spin configurations. By 
analogy with proton ordering in ice IhEJa, the spin con- 
figuration with two spins pointing inward and two spins 
pointing outward for each tetraheron corresponds to the 
'ice-rule' for these compounds. Accordingly, Ho2Ti207 
and related compounds such as Dy2Ti207, which also 
possesses FM Ising spins on a pyrochlore lattice, were 
named 'spin ices'. 

The first naive theoretical model of the spin ice materi- 
als placed Ising spins on a pyrochlore lattice, constrained 
to point along the local (111) direction- with only nearest- 
neighbor FM exchange interactionsBB. However, it was 



soon deduced that the nearest-neighbor exchange inter^. 
action in these compounds is actually antiferromagneticu 
and that the effective nearest-neighbor interaction be- 
comes ferromagnetic due to the long-range dipolar inter- 
action arising from the large moments of the Ho 3+ and 
Dy 3+ ionsQ. The physical model which correctly incor- 
porated all of these interactions is called the dipolar spin 
ice model. Physical quantities (neutron scattering cross 
sections and specific heat curves) calculated from Monte 
Carlo simulations of this model are in excellent agree- 
ment with the corresponding experimental dataLTu. 

Most experimental studies on the spin ices ■La.jdaie 
have been performed on polycrystalline samplcsomj 
Among these various experiments, the most dramatic and 
direct evidence for realization of the spin ice state is the 
residual entropy equal to Pauling's entropy: ^i?ln| = 
1.68 J/K 2 mo]o. Ramirez et al. obtained an entropy very 
close to this value (~ 1.9(2) J/K 2 mol-Dy) by integrat- 
ing the experimental ekctronic specific heat C e /T of 
Dy2Ti207 above 0.25 KB. This residual entropy can be 
removed by applying an external magnetic field, since 
the magnetic field can lift the macroscopic degeneracy 
by stabilizing a specific spin configuration with lower 
magnetic energy (Fig. 1), as proven by neutron diffrac- 
tion measurement on single crystals of Ho2Ti207 and 
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FIG. 1. Spin configurations with the magnetic field along 
(a) the [100], (b) [110] and (c) [111] axes. For the [111] direc- 
tion the figure depicts the state in which the ice-rule breaks, 
(d) Three dimensional axes in order to give three dimensional 
perspective of the direction in which the field is applied. 

Moreover, the electronic specific heat, C e (T)/T, of 
Dy2Ti207 polycrystals exhibit three peaks at 0.34 K, 
0.47 K and 1.12 K independent of the magnitude of the 
magnetic fieldO. Due to the freedom of the decoupled 
spins of Fig. 1(b), it is speculated by Ramirez et al. that 
these peaks are attributable to crystal grains with their 
[110] axis aligned parallel to the field. 

The above experimental facts suggest that the spin 
ice compounds have significant magnetic anisotopy, and 
therefore one must control the field direction carefully 
when studying their physical properties. In particular, 
the competition between the applied magnetic field and 
the spin-spin interactions that drive the ice-rule obeying 
degenerate ground state is significanty affected by the 
crystal orientation in the field . However, to date there 
has been no experimental report detailing the behav- 
ior of this anisotropy using single crystals of Dy2Ti2C>7. 
In this paper, we report such measurements and com- 
pare our results with the results predicted by the dipo- 
lar spin ice model. We show experimental magnetiza- 
tion results at 1.8 K along three characteristic axes of 
this compound, namely the [100], [110] and [111] axes, 
and the corresponding numerical results calculated by 
Monte Carlo simulations. As we will discuss, agreement 
between the two is excellent. In contrast, a similar set of 
measurements carried out by Cornelius and Gardner on 
Ho2Ti 2 07 failed to find the expected magnetic anisotropy 
associated with this materials. 

In addition, it has recently been theoretically predicted 
that for the dipolar spin ice model and, consequently, for 
DV2TI2O7 and Ho2Ti207, a first-order phase transition 
should occur at around 180,-lqK (= T c ) which fully re- 
leases Pauling's ice entropyfij. Melko et al. predicted 
the transition by considering a spin loop move instead 



of single spin flips within their Monte Carlo simulation, 
as illustrated in Fig. 1 of Ref. [l5]. However, few ex- 
plicit experimental searches for this transition have been 
reportedMES. Hence, the second motivation for this pa- 
per is to investigate single crystals of Dy2Ti207 using 
susceptibility measurements down to 60 mK, in order to 
search for the predicted magnetic ordering. 

It is hoped that this experimental and theoretical joint 
study on Dy2Ti2 07 can help answer some current issues 
of debate, and help develop deeper understanding of the 
nature of the spin ice materials. 

II. EXPERIMENTS AND CALCULATIONS 

Single crystals of Dy 2 Ti 2 07 were grown by the float- 
ing zone method with an-infrared furnace equipped with 
double-elliptical mirrors£3. Before the single crystal 
growth we prepared a polycrystallinc feed rod by thor- 
oughly mixing Dy 2 3 (99.99%) and Ti0 2 (99.99%) pow- 
ders with the molar ratio of 1 to 2 and sintering the mix- 
ture in an alumina crucible at 1000°C for 100 hours in 
air. The typical growth conditions are 29.0 mm/h for the 
feed speed and 30.0 mm/h for the crystal growth speed. 
To avoid oxygen deficiency, we grew the single crystals 
in 0.3 MPa of O2 atmosphere. In addition, we annealed 
the as-grown crystals at 800°C for 100 hours in O2 gas 
flow. The color of the annealed crystals is homogeneously 
translucent yellow, while that of the as-grown crystals is 
lighter yellow with some inhomogeneous color distribu- 
tion. The sizes of samples used for magnetization and ac 
susceptibility measurements are 0.7x 0.7x 2.0 mm 3 and 
3.0 x 1.0 x 0.2 mm 3 , respectively. By performing powder 
x-ray diffraction measurements on the annealed crystals, 
we found that the scattering pattern was consistent with 
a cubic unit cell with a = 10.124(5) A without any ex- 
tra peaks. We determined the principal axes using Laue 
pictures. 

We measured magnetization between 1.8 aad 40 K for 
fields up to 5 T with a SQUID magnetometertS. We mea- 
sured ac susceptibility by a mutual-inductaace method 
down to 60 mK using a dilution refrigeratorlij. 

For magnetization measurements, we applied the mag- 
netic field in the plane perpendicular to the longest side 
of the crystal, because all of the principal axes are con- 
tained within this plane. The fields-misalignment is esti- 
mated to be within a few degreeai3. Thus, we estimate 
the demagnetization factor n as 0.4 for this crystal. The 
demagnetization effect can be neglected for ac suscepti- 
bility measurements, since the ac field was applied par- 
allel to the longest side, which is parallel to the [100] 
crystal axis. 

In addition, we performed Monte Carlo simulations for 
magnetization at finite temperatures on the dipolar spin 
ice model. The detailed calculation techniques are de- 
scribed in Ref. (?]. We adopted the following Hamilto- 
nian: 
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Here is the Ising moment vector of magnitude of 
unity at site i pointing along the local Ising axis Zj. 
The quantity JT ff is the effective internal magnetic field 
which directly couples to each spin. For the Hamiltonian 
coefficients, we used —1.24 K for the nearest-neighbor ex- 
change interaction J nn = J/3 and 2.35 K for the nearest- 
neighbor dipolar interaction D nn = 5D/3. We took into 
account the long-range nature of the dipolar interactions 
through the Ewald summation technique. A standard 
Metropolis algorithm was used with a conventional cubic 
unit cell for the pyrochlore lattice containing 16 spins. 
Simulations of the specific heat obtained using these val- 
ues are in very good agreement with the experimental 
specific heat reported for Dy2Ti2C>7. The simulation size 
was L = 3 (432 spins) . When calculating the magnetiza- 
tion curves, a number of different simulation sizes were 
tried, and little variation resulted. 



III. RESULTS AND DISCUSSION 
A. Magnetic Anisotropy 

In Fig. 2 we summarize the magnetization curves at 
1.8 K along the three characteristic axes. We took the 
horizontal axis as the effective magnetic field, /logoff = 
Mo-ffcxt — nM exp . We obtained nearly the same data at 2 
K (not shown). We note that no hysteresis is observed 
between the upward and downward field sweeps. 

Saturated moments along the [100] and [110] axes agree 
with the expected values for Ising anisotropy along the 
local (111) direction and for effective FM coupling be- 
tween the nearest-neighbor Dy ions. The value for the 
[100] axis is gjJ(l/y3) — 5.77 ^e/Dy (two spins point 
in and two spins point out: Fig. 1(a)) and that for the 
[110] axis is g j J {{^J2jl x 2) / A} = 4.08 Me/Dy (one spin 
points in, one spin points out and two spins are free: 
Fig. 1(b)). Here, gj represents Lande's g factor, J cor- 
responds to the total angular moment of Dy 3+ ion, and 
Hb represents the Bohr magneton. For the free Dy 3+ ion, 
gjj is expected to be 10 /Us/Dy. Slight deviation from 
the expected values may be due to misalignment effect 
as we discuss below. 

Below 0.5 T, the magnetization along the [111] axis 
already reaches the value gjJ{(l + 1/3 x l)/4} = 3.33 
Mb /Dy expected for the saturated state in which the ice- 
rule completely holds. In this state, three of the spins 
per tetrahedra are coupled to minimize their energy with 
the applied magnetic field along the [111] direction, but 
one spin chooses to oppose the field and minimize its 



energy with respect to the effective nearest-neighbor in- 
teraction, remaining in the ice-rule. For quite large fields 
(H c g » J nn + D nn ), this configuration begins to cost 
a large amount of energy, and the field-opposed spin 
chooses to flip, minimizing its energy with respect to the 
applied magnetic field, but breaking the ice-rule. Ac- 
cordingly, the magnetization agrees with the value ex- 
pected for the three-in one-out spin configuration (Fig. 
1(c)), gjJ{(l + 1/3 x 3)/4} = 5.00 /x B /Dy. We infer 
that the gain of Zeemann energy exceeds that of the 
effective nearest-neighbor FM interaction above 0.5 T: 
(l/3)gjJfiB(fJ'oHes) = 1.12 K, in good accordance with 

J e g = J nn + Am = l.llK. 
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FIG. 2. Magnetization curves at 1.8 K along the [100], 
[110] and [111] axes. Solid lines correspond to the expected 
saturated moments for the two-in two-out spin configuration 
for all the direction: dotted line corresponds to that for the 
three-in one-out spin configuration. 

For fields up to 0.5 T, the moments along all the axes 
nearly reach the expected values for the corresponding 
spin-ice configurations in a magnetic field. This exper- 
imental fact indicates that the magnetic anisotropy in 
Dy2Ti207 is indeed compatible with the spin-ice model. 
In addition, the breaking of the ice-rule observed along 
the [111] axis confirms that the strong Ising anisotropy 
along the local (111) direction dominates in this com- 
pound. 

A corresponding field dependence of magnetizations is 
reported by, Cornelius and Gardner for single crystals of 
Ho2Ti2070. However, their data above 2 T is totally 
unexpected, because the magnetizations along all the di- 
rections reach approximately 5.9 /zb/Ho. Their behavior 
cannot be explained solely by a crystal misalignment ef- 
fect, since the effect of misalignment eaasts in our exper- 
iment but results in different bchavioiEa. 
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In Fig. 3 we show magnetizations at 1 and 2 K calcu- 
lated by the Monte Carlo simulations. Each solid sym- 
bol represents the magnetization along the correspond- 
ing held direction of the experimental results. The open 
squares are for the magnetization at 2 K along the [110] 
axis with intentional 0.5 degree misalignment toward the 
[100] axis. The open diamonds are for the magnetization 
at 1 K along the [111] axis. 

Clearly the Monte Carlo simulations at 2 K reproduce 
our corresponding experimental magnetization curves 
very well, though the absolute values do not completely 
coincide with each othere3. This excellent agreement be- 
tween experiments and calculations implies that one may 
treat this compound as an ideal model system for the 
study of geometrical frustration in the pyrochlore lattice. 
Moreover, it shows that the assumption of an Ising dou- 
blet with leading wavefunctions \ip-\-) an d |^— ) with \mj) 
decompositions \ip±) ~ \mj = ±15/2) is a very good 
approximation, in goocL-agreement with crystal-held cal- 
culations for Dy2Ti207E3 adapted from the crystal-held 
parameters extracted for Ho2Ti 2 OTO. 



deviation from the expected value, 4.08 /le/Dy, can be 
explained by assuming a crystal misalignment within 0.5 
degree from the [110] axis. 

The derivative of the experimental and calculated mag- 
netizations at 1.8 K with H e g along the [111] axis are 
shown in Fig. 4. The derivative of the experimental 
magnetization, dM cxp /dH c g, exhibits a definite kink at 
around 0.5 T and quite a broad peak at around 1 T. 
This kink signals the crossover from two- in, two-out spin- 
ice conhguration to three-in, one-out conhguration. The 
derivative of the calculated magnetization, dM cgu \/dH c g, 
also exhibits a local minimum at around 0.5 T and a 
broad peak at around 1 T. In addition, the behavior be- 
low 0.5 T nearly coincides with the experimental results. 
This good agreement between experiments and calcula- 
tions indicates that the behavior of the dipolar spin ice 
model in the spin-interacting regime (i. e. at low helds) 
behaves almost the same as the real materials in the same 
regime. The difference between them is most likely due 
to the experimental misalignment effect, since the differ- 
ence became smaller when we measured another single 
crystal by adjusting the alignment more accuratelyEJ. 
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FIG. 3. Magnetization curves calculated by Monte Carlo 
simulations of the dipolar spin ice model. The simulation 
size was L = 3 (432 spins), and parameters were chosen for 
DyzTiaOr. 

The calculated magnetization with the held along the 
[110] axis with the intentional 0.5 degree misalignment 
deviates substantially from the aligned value at higher 
helds. This is due to the acute angular sensitivity specific 
to held application along the [110] axis, since slight devi- 
ation of the held from the [110] axis can hx the spin direc- 
tions of the 2 free spins normal to the [110] axis in every 
tetrahedron (Fig. 1(b)). The increase in magnetization is 
attributable to the held-coupled component of these spins 
still satisfying the ice-rule. Therefore, the experimental 




FIG. 4. Derivative of the experimental and calculated 
magnetizations at 1.8 K with iiqH c s along the [111] axis. 

The calculated broad peak becomes sharper at 1 K and 
correspondingly a clear plateau between 0.25 and 0.45 T 
emerges in the magnetization curve at 1 K. The plateau 
definitely suggests that an ice-rule obeying spin configu- 
ration is more stable than the three-in one-out spin con- 
hguration at lower temperatures and at lower helds. Such 
plateau is-sdso reported in Ho2Ti207 at 1.8 K between 0.4 
and 1.6 Tli-a. In Ho2Ti 2 07, the appearance of the plateau 
clearly corresponds to the temperature of the maximum 
in the specific heat. n This characteristic temperature is 
1.97 K in Ho2Ti20va and corresponds to the onset of 
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strong spin-ice correlations in the material. The corre- 
sponding characteristic temperature in Dj^T^Oy is 1.2 
KB. The results that a clear plateau appears at 1 K in 
our calculations but not above 1.8 K in our experiments 
as well as our calculations are reasonable compared with 
the characteristic temperature for Dy2Ti20y. 

In Fig. 5 we show the temperature dependence of 
reciprocal dc suscpeptibility for various magnetic fields 
along the [100] axis. The solid line denotes the Curie- 
Weiss fitting between 16 and 40 K for the data at 0.01 
T: 
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(3.1) 



where x = M/H is magnetic susceptibility, Na is the 
Avogadro number, p e g is effective Bohr magneton, and 
/cb is the Boltzmann constant. Again, H has been cor- 
rected for the demagnetization factor. 
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FIG. 5. Temperature dependence of reciprocal dc suscpep- 
tibility for various magnitudes of magnetic field along the 
[100] axis. 

For a given temperature, the ratio H / M increases with 
an increasing magnetic field. This is because the magne- 
tization at higher fields saturates compared with that at 
lower field especially at lower temperatures. 

We obtained positive Curie- Weiss temperatures, 
1.16(5), 1.31(5) and 1.28(5) K from the data at 0.01 T 
along the [100], [110] and [111] axes, respectively. This 
is consistent with the existence of the effective nearest- 
neighbor FM interaction, mainly originating from large 
dipolar interaction between Dy ions, and consistent with 
the previous results using polycrystalline samplesB'£3. 
Moreover, this provides further validity in using the value 
for the effective FM interaction used in our calculations, 

Jeff = J nn + Am = 1-11 Kc3. 



We also obtained the effective Bohr magneton 10.2(1), 
10.0(1) and 10.0(1) MB/Dy for the [100], [110] and 
[111] axes, respectively, from the Curie- Weiss fitting. 
These values are nearly equall to the value p c ff = 
gjyJ[J + 1) = 10.6 ^e/Dy expected for the completely 
localized spin of Dy 3+ ion. 



B. Ac Susceptibility and Search for the 180 mK 
Transition 

Considering the above excellent agreement between 
our experiments and calculations, we may expect the re- 
alization of the theoretically predicted first-order phase 
transition in the real materialtJ. In Fig. 6, we show the 
ac susceptibility Xac down to 60 mK. In order to keep the 
sample in zero dc field, we surrounded the sample coils 
with high-permeability paramagnetic tube. An ac field 
of 0.05 mT was applied along the [100] axis. This am- 
plitude satisfies the condition, H cxt / D nn ~2x 10~ 4 <C 



0.08 ~ T c /D 

nni suggested in Ref. |l 5j . We measured ^ ac 
at 1 kHz on cooling and at 17 Hz on warming after keep- 
ing the crystal at 60 mK for 1 hour. We normalized the 
data by dividing the raw data by the measuring frequecy. 
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Ac susceptibility measured down to 60 mK. 



The in-phase component x' measured at 17 Hz exhibits 
a broad peak at around 1.2 K and steep decrease below 
this temperature. In addition, it diminishes below about 
0.5 K. x' at 1 kHz and the out-of-phase, dissipasive com- 
ponents x" a f both frequencies also exibit the similar 
temerature dependence. These data are consistent with 
those obtained from polycrystalline Dy2Ti20vll3. 

The absence of any change in the susceptibilities below 
0.5 K indicates that the spins on each tetrahedron in the 
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system become firmly locked in a random ice-rule config- 
uration at about 0.5 K, below which the specific heat also 
diminishes as shown in the previous study by Ramirez et 
alu. Further, our data suggest that the predicted first- 
order phase transition does not occur, at least within the 
limited experimental measurement time scale. Because 
the 'freeze in' temperature is well above the predicted 
ordering temperature of 0.18 K, it is plausible that there 
is insufficient local spin dynamics left in the system that 
can promote the ordering associated with the nonlocal 
dynamics of spin loop move in the Monte Carlo simula- 
tions of Rcf. [l5| We therefore speculate that the sys- 
tem has essentially completely fallen out of equilibrium 
at around 0.5 K; below this temperature one does not ex- 
pect to recover equilibrium in a time that is reasonably 
accessible by experiment. 

The temperatures of the maxima in \' an d x" increase 
with increasing frequencies. The dependence is ascrib- 
able to long magnetic relaxation-time of spin lock-in for- 
mation as Matsuhira_et al. discussed in their previous 
polycrystalline studylifj. This also supports our interpre- 
tation of the absence of ordering down to 60 mK. In fur- 
ther studies, it may be possible that slight dilution of 
Dy with non-magnetic elements such as Y may bring the 
static ordering, in analogy to the transition in water ice 
from the Ih phase to the long-rapge ordered crystal Ixi 
phase via substitution with KOHc3. 



IV. CONCLUSION 

In conclusion, we reported magnetization and ac sus- 
ceptibility of Dy2Ti2C>7 using single crystals. The clear 
magnetic anisotropy indicates the strong Ising anisotropy 
and the effective FM interaction between the nearest- 
neighbor Dy ions, as well as the dictating ice-rule at 
low temperatures. These effects have been quantita- 
tively confirmed by both experimental and theoretical 
approaches in the present study. 

No indication of the predicted long-range magnetic or- 
dering was observed. This may be due to differences in 
time scale between the dynamics of the real system and 
the simulation. Futhcr investigations, such as specific 
heat measurements, or dilution with non-magnetic im- 
purities may clarify the attainability of the real ground 
state of Dy 2 Ti 2 7 . 
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